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Abstract. We report the first measurements of the inelastic spin exchange collision 
rate between the Hyperfine (HF) levels of antiprotonic helium (pHe + ). We measure 
the time dependent evolution of the (37, 35) substates to obtain an inelastic collision 
rate which qualitatively agrees with recent theoretical calculations. We evaluate these 
results by using the obtained rate as a parameter in a rigorous simulation which we 
then compare to to previously measured data. We find that our measurement slightly 
underestimates the collision rate and therefore conclude that the actual value most 
probably falls within the upper, rather than lower, limit of the error. 
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Figure 1: (a) Level diagram of p 4 He + where the arrows indicate the radiative cascade 
towards the nucleus, (b) Hyperfine splitting of an (n, /) state of p 4 He + . The wavy lines 
denote allowed Ml transitions that can be induced by an oscillating magnetic field. 
From [1]. 

1. Introduction 

Antiprotonic helium pHe + is formed when an antiproton p interacts with a helium atom 
at or below the ionization energy (~ 25 eV) [2]. The p can become captured so that 
it precesses around the helium nucleus He ++ [1,3,4]. When this happens, one of the 
electrons e~ is ejected. Because of its mass, the p is most likely to occupy an orbit with 
principle quantum number n = n = M* /m c ~ 38 [5], where M* is the reduced mass 
of the antiproton-helium nucleus system and m e is the electron mass. It precesses in a 
semi-classical orbit while the electron remains in a Is quantum mechanical cloud. 

Because of their overlap with the nucleus, the majority of captured antiprotons 
annihilate within picoseconds with one of the nucleons in the He ++ nucleus [4] . However, 
approximately 3% occupy metastable states, so called circular states in the region of 
n = 32-40 and vibrational quantum number v = 0-3 (where v = n — I — 1 and I is the 
angular momentum quantum number). Since the neutral system retains one electron, it 
is protected from external atoms by the Pauli exclusion principle [4]. Additionally, the 
presence of the electron removes the / degeneracy for the same n, therefore protecting 
it against Stark mixing. The Auger decay of the remaining electron is suppressed by 
the large ionization energy compared to the n — > n — 1 level spacing of ~ 2 eV. Thus 
only one decay channel remains and the antiprotons in these states undergo a radiative 
cascade through (n, I) — > (n — 1, I — 1) states, each with lifetimes in the order ~ 1.5 /xs, 
see figure da! 

A hyperfine (HF) splitting [4], caused by the interaction of the e~ spin S c with the 
p orbital angular momentum L, results in a doublet structure of the order z^hf = 10- 
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Figure 2: (a) Schematic view of the primary level splitting of pHe for the unfavoured 
electric dipole transitions. The state drawn on the right is the radiative decay dominated 
parent (n, L), and state on the left is the Auger decay dominated daughter (n', L'). The 
laser transitions, from the parent to the daughter doublets, are indicated by the arrows 
/ + and /~. (b) Laser resonance profile of the (n,l) = (37, 35) to (38, 34) transition 
showing the / + and /~ transitions, from [7] 



15 GHz. A further splitting of each HF state results in a superhyperfine (SHF) structure 
(^shf = 150-300 MHz), caused by the interaction of the p spin Sp with F = L+S c . There 
exists therefore a quadruplet substructure for each (n, I) state as shown in figure [TBI 
The theoretical framework for the level splitting has been developed by Bakalov and 
Korobov [6]. 

The HF doublet is described by the quantum number F = L + S e with components 
F + = L + ~ and F~ = L — |. The SHF quadruplet is described by J = F + Sp with 
components J~ + = F~ + ~, J~~ = F~ - §, J ++ = F + + \ and J + ~ = F + - |. Between 
these sub states an electron spin flip can be induced by two Ml transitions and u^ F 
(see figure [Tbl) : 



F + + ^ = L + 1^ J-+ = F~ + ^=L, (la) 

F + -- = L <-> J — F~ — - — L — 1. (lb) 
2 2 v ; 

Electric dipole transitions (El) between different levels of the cascade can be 
induced with resonant laser light [8—11]. There are two types: favoured, Av = 0; 
(n, I) — > (n — 1, 1 — 1), and unfavoured, Av = 2; (n, I) — > (n + 1, I — 1) [4]. The dipole 
moment of the unfavoured transitions is an order of magnitude smaller than that of the 
favoured. To the first order, atoms occupying the F~ doublet of the radiative decay 
dominated state are transferred to the F'~ doublet of the Auger dominated state, while 
those occupying F + are transferred to F' + , shown in figure [2H These transitions are 
labelled / + and /~ respectively and the difference between them Af. The unfavoured 
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transitions have Af = 1.5 — 1.8 GHz, while the favoured have Af < 0.5 GHz. The 
Doppler broadening at the target temperature is ~ 0.3 — 0.5 GHz, therefore only the 
unfavoured HF laser transitions can be well resolved, see figure l2bl 

2. Motivation and Method 

A precise measurement of the pHe + HF splitting [7, 12, 13] is of great importance for 
rigorously testing three-body quantum electrodynamic (QED) calculations [6, 14-16], 
leading to a determination of the antiproton spin magnetic moment and a test of CPT 
invariance [17]. For an understanding of the collision processes between pHe + atoms 
and the He atoms of the medium, a comparison between experiment and theory can be 
equally useful. Both the elastic T e and inelastic Tj collision rates can have significant 
systematic effects on experimental results. Elastic collisions contribute to a shift and 
a broadening A — iT e , while inelastic collisions result in a spin exchange, therefore a 
change of state. 

A clear understanding of the collision processes was essential for the interpretation 
of the El spectral lines [18], a similar study has been performed for the Ml transitions. 
To measure A — iT e , microwave resonant profiles were scanned at various target gas 
densities. The line width is limited by the Fourier transform of the microwave pulse 
length but observation of a larger width would be evidence of a collisional broadening. 
Likewise a density dependent change in the transition frequencies would be evidence of 
a collisional shift. 

The results of elastic collisional studies have been presented in previous 
publications [7, 19] and indicate that r e is small because the dominating broadening 
effect is found to be from the Fourier transform of the microwave pulse length. Korenman 
predicts that T e ~ 2.5Fj [20] which means that if T e is smaller than first predicted then 
so must Tj. 

The inelastic collision rate was determined by measuring the time dependence of 
the F + population. Two narrow-band lasers were tuned to the / + transition between 
the radiative decay dominated parent state (n, L) = (37, 35) and the Auger decay 
dominated daughter state (38, 34) shown in figure [2j The second was delayed by a time 
T = 50 — 2000 ns from the first. 

The p annihilation products were detected with Cherenkov counters as a function of 
time. The metastable tail, where the radiative decay dominated states cascade towards 
the nucleus, was recorded as background. Because of its short lifetime (~ 10 ns), the 
laser resonant transfer to an Auger dominated decay state results in a sharp peak in 
annihilations events which stands out against the background, shown in figure [3l The 
ratio between this peak area to the area under the entire spectrum (peak-to-total) is 
proportional to the population transferred with the laser. The peak-to-total of the first 
and second laser annihilation peaks are represented by r\ and r2, respectively. 

The experiment was performed in two different modes: 1) where both lasers were 
fired, f + -f + and 2) where only the second laser was fired, 0-/ + , both of which are shown 
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Figure 3: Laser stimulated annihilation peaks against the exponential decaying 
background of the other metastable states' populations. Mode 1): f + -f + when both 
lasers are fired and mode 2): 0-/ + when only the second laser is fired. The peak in 
Mode 2) is larger than the second peak in Mode 1) because no previous depopulation 
of the state has been induced. 



in figure [HI Mode 1) contained all the information about Tj while mode 2) was required 
to extract information about the refilling from higher states, also contained in 1). 

The method employed in [13], [7] and [17] to determine the HF splitting included 
a microwave pulse between the first and second lasers of mode 1). By scanning the 
microwave over a range of frequencies and measuring the dependence of r%, the z/^p 
and Vft F resonances were found. The maximum achievable signal for a range of laser 
delays was determined by fixing the microwave frequency to one transition, say u^ F , 
and monitoring r 2 while scanning the microwave power. 

3. Apparatus 

The experiment was performed at CERN's Antiproton Decelerator (AD), which 
delivered a pulse of 1-4 x 10 7 antiprotons with a length of 200 ns (FWHM) and an 
energy E = 5.3 MeV at ~ 90 s intervals. Antiprotonic helium was formed by stopping 
antiprotons in a gas target at a temperature of 6.1 K and a pressure p = 150-500 mbar 
(number density 1.7-6.2 x 10 20 cm" 3 ). 

Charged pions were produced by antiproton annihilations in the helium nucleus 
and could be detected by Cherenkov counters. The signal was amplified by fine-mesh 
photomultipliers (PMTs) and the resulting analog delayed annihilation time spectrum 
(AD ATS) was recorded in a digital oscilloscope (DSO). The PMTs were gated off during 
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the p pulse arrival so that only the 3% metastable tail was recorded [21]. 

Two pulse-amplified laser beams were produced by splitting a continuous wave 
(cw) laser beam of wave-length 726.1 nm into two seed beams [7]. These were pulsed 
by amplifying the seeds using dye filled Bethune cells pumped by two pulsed Nd:Yag 
lasers, the second delayed by time T after the first. The pump beams were stretched 
so that the two pulse lengths were of the order ~ 15 ns [11] and therefore similar to 
the Auger decay rate ensuring a high depopulation efficiency The maximum emitted 
energy fluence at the target was ~ 30 mJ/cm 2 with a spot diameter of 5 mm. 

To measure the HF transitions, a microwave pulse was produced by a vector network 
analyzer (Anritsu 37225B) referenced to a 10 MHz GPS (HP 58503B) satellite signal 
and amplified by a pulsed travelling wave tube amplifier (TMD PTC6358). A cylindrical 
resonant microwave cavity with central frequency uq = 12.91 GHz provided the desired 
shape for the field (TMno mode) at the target. To cover the Av ~ 100 MHz microwave 
scanning range, the cavity was over-coupled to the wave guide so that its loaded quality 
factor Ql was ~ 100, where Av = /o/Ql [22]. Most of the power was reflected back 
towards the amplifier and absorbed by an isolator. An antenna was connected to the 
cavity to monitor the field so that the desired power could be achieved by controlling 
the amplification of the pulse. 



4. Analysis 

4.1. Mode 1) (f+-f+) 

When both lasers were fired, the second delayed by a period T after the first, the 
normalized peak-to-total ra/ri was plotted as a function of T. The data were fitted 
with a function derived from the integral of the following two equations: 

dp H 



dt 
dp_ 



9+(t) - (A+_ + 7r)p+ + A_+p_, (2a) 

g-(t) - (A_ + + 7r )p_ + X+-P+, (2b) 

where p± is the population density of the HF states, and g±[t) describes the refilling 
rate from the higher lying states. The relaxation rates from p + to p_ and p_ 
to p + are represented by A + „ and A_ + respectively. The radiative decay rate is 
7r = 7.149 x 10 5 s" 1 [4], see figure H 

If t — is the time when the first laser is fired then the relative population of the 
two levels at t < is 

P± ~ 2(2L + 1)' (3) 

where F + = L + | and F~ = L — |. The signal from the first and second laser are 
therefore 

n = Io^-jet, (4a) 
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Figure 4: Energy level diagram of part of the cascade showing the refilling g±(t) from 
above states, decay to lower states 7 r and the relaxation collision rate A±^ T . 



r 2 (t) = I p+(t)e 2 , {4b) 

where Iq is a normalization factor and e\ and e 2 are the laser depopulation efficiencies 
for the first and second laser. The overlap of the laser peaks, whereby the F~ transition 
is partially induced due to the Doppler broadening when the laser is tuned to the f + 
transition, is considered negligible in this analysis. 

At t = and since L ^> 1, the initial populations are 

P+(0) = ^1(1 - £l ) « i(l - (5a) 

"- (0 > = 2ZTT 4 < 54) 

It can also be assumed that A + _ = A_ + = Tj and g + = g_ = g. So (12 a\) and (j26|) can 
be written as follows 

P+ (t) = ie-**(l - | + |e~ 2r ^) + 1 fe-^y f h(t')dt', (6a) 

J 

P-(t) = - | - |e- 2r ^) + \ te-^vMt'W, (66) 

where g = fifh(t), of which /if is a constant associated with the rate of filling from above 
states and h(t) is the filling function normalized by h(0) = 1. The simplest assumption 
for h(t) is h(t) = e _/i °* where /j,q is a decay rate associated with the population of the 
above states. This simplification was necessary to achieve an unambiguous result due 
to the limited amount of data. So the solution to the integral in (RJqj) and fl66|) becomes 

F(t) = N f e-^-t')h{t')dt' = /i/(e " 7rt ~ e ""° f) . (7) 

J0 PO ~ lr 
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At T < 1000 ns the filling can be assumed constant llq = 0, and thus 
F(t) 



/*/(l-e->*) 



Substituting (jSJ) into (|6aj) . then flo r aj) into (j4Dl) and normalizing over (|4qp gives 



r 2 (*) 

n(t) " 

Mode (o-/+j 



£2 
El 



2 2 



lr 



(8) 



(9) 



Only firing the second laser, the population decay and refilling of the state can be 
measured. In this regime r2 was normalized to the average r± from mode 1). This was 
performed so that both sets of data could be plotted on the same scale and compared 
adjacently. There are no inelastic collision terms because no asymmetry is created, 

r 2 (t) 

" 1 1 ~) ^ ~ 



£2 



(10) 



4-3. Numerical simulation 



The (37, 35) — > (38, 34) laser transition was numerically simulated by evolving the 
optical Bloch equations, obtaining a maximum depopulation efficiency of 70%. 

The microwave transitions between the HF substates were determined by evolving 
(II lap , derived from the optical Bloch equations, which can be written as two independent 
4x4 matrices to handle the z/^ F and transitions separately. However, collision 
induced p spin flips result in v$ UF and z/^" HF transitions. Thus the population evolutions 
of the J h and J ++ states become dependent on those of the J and J~ + states. The 
resulting simultaneous equation has the solution of the form of an 8 x 8 matrix. An 
additional two dimensions were added to simulate the refilling from above states: 

d 



dt' 



= Mp 

{ P-+\ 
P++ 

Px+ 

Py+ 
P— 

P+- 

Px- 

Py- 

V Pu39 j 



[lla) 



(lib) 
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where p_ + ,p ++ ,p , and p + _ represent the different time dependent populations of 

the four SHF states of the (n, I) = (37, 35) state. In a field free region Q m = 0, these 
populations simply decay radiatively at a rate of 7 r = 7.149 x 10 5 s _1 to the (36, 34) state 
and 0.0086 x 10 5 s _1 to the (37, 34) state [4]. When there is a population asymmetry 
and an external oscillating magnetic field is present, transfer between the states can be 
observed. The complex dependency of the transitions is represented by p x ± and p y ± for 
the real and imaginary parts respectively. Broadening effects are dependent on both the 
radiative decay rate 7 r and the elastic collisional frequency T e : = lr + r e . 

The four SHF states are refilled as the upper states, (38,36) p U 3$ and (39,37) p„39, 
decay at a rate of 7 u38 = 6.55 x 10 5 s _1 and 7 u39 = 5.88 x 10 5 s" 1 into the lower 
(37, 35) state, as part of the cascade [4]. The initial populations of these states have 
been experimentally measured [18,23]. Through inelastic relaxation collisions Tj the 
atoms return to an equilibrium, the variable 7 C is defined as 7 C = 21^ + 7 r . Collisions 
which result in the spin flip of more than one particle are ignored [20]. 

The offset between the microwave frequency u M and transition frequencies z/hf is 
represented by Alu = 2tx(vm — ^hf±)- The Rabi frequency is dependent upon the 
magnetic field strength B and the atom's magnetic dipole moment fi m 

p rn B(x,y,t) 



h 



p m = (n', L', F', J', m\p M \n, L, F, J, m), 
which can be calculated using Wigner's 3-j and Racah's 6-j coefficients 



(12) 
(13) 



(n,L,F , J', m\p\n, L, F + , J, m) = (— 1 



v J'+m 




3-j 



x v/(2J+l)(2J' + l)<! F ~ J p+ \ 



(14) 



xv/(2F+ + l)(2F- + l)<i j + F ~ \ 



6-j 



9eM-\s e \-) 



where n, L, F, J, m are the corresponding quantum numbers. 
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Figure 5: Simulated microwave resonance profile, where T = 350 ns Tj = 3.4 x 10 5 s 1 , 
fitted with the same function as the experimental measurements presented in [7] and [17] 



Equation filial) was solved over a range of 60 equally spaced frequencies vm 
to simulate a microwave resonance profile measurement. The population positional 
distribution and magnetic field variance were modelled with a Monte Carlo positional 
simulation. The magnetic field distribution at the target is dependent on the cavity and 
varies spatially with respect to the radial r and angular cylindrical <fi co-ordinates [22]. 
Apart from edge effects the cylindrical component z is constant. The radial B r and 
angular components of the magnetic field are given by 

B r {r,(t>) = 5o^^sin(0), (15a) 

B 4) {r,<l>) = B Q J' 1 {kr)co&{<j>), (155) 

where k is the wave number and J\ is the Bessel function of the first kind. The stopping 
distribution of the pHe + is assumed to be Gaussian in both the z and r planes [22] . 

5. Results 

Data were measured at three different target pressures p = 150 mbar, 250 mbar and 
500 mbar resulting in a total of six graphs; three for mode 1) and three for mode 
2). These were plotted side by side and fitted simultaneously with (J9j) and (JTOl . The 
variables fif, t\ and Tj were common for all pressures, where the latter was weighted 
proportionally to the target gas density and €2 was left free for different target densities. 
Other, more complex fit functions were also attempted. These varied to include the 
population evolution of the upper levels and left e\ free for different target densities. 
However, the introduction of more parameters limited convergence and put emphasis 
on the refilling processes. It was found that the simplest function provided the most 
sensitivity to IV 
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Figure 6: Simultaneous fitting of the data for p = 150 mbar, 250 mbar and 
500 mbar. Mode 1) (f + -f + ) was fitted with @ and mode 2) (0-/+) was fitted 
with (flOl) . The refilling rate /i/ was a common parameter for all target densities and the 
laser depopulation efficiencies t\ and €i were common parameters for different modes 
measured at the same density. 
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Figure 7: Dependence of the collision induced relaxation rate Tj on the helium target 
density. 



A graph of the data fitted with ([9]) and ( flTJl) is shown in figure EJ The collision 
induced relaxation rate has been determined from the fitting parameters and plotted in 
figured The numerical values are displayed in table [1] except for \ij and e\ which were 
determined to be (5.2 ± 0.2) x 10 5 s _1 and 55% respectively. 

The errors associated with these relaxation rates have been inflated by the square 
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Figure 8: The optimum signal-to-noise ratio for each time delay T at p = 250 mbar. The 
circles (•) represent the experimental results from Pask et. al. [7] while the triangles 



(a) represent the simulated data. The larger signal results from Tj 
the smaller from T, = 4.5 x 10 5 s _1 



2.3 x 10 5 s" 1 and 



root of the reduced chi squared Xred ~ 10 of the fit but still remain 2-3 sigma less 
than the most recent theoretical calculations which predict Tj = 6.2 x 10 5 s -1 [24] for 
p = 250 mbar. The laser depopulation efficiency is also revealed to be smaller than the 
70% predicted in Section 14731 The fluctuations of €2, that can vary by as much as 15%, 
are put down to the fact that data of different densities were measured on different days. 
The fluctuations are therefore most probably caused by changes in the overlap of the 
two lasers with the p beam and each other. 

To be certain that the fit provided a realistic determination of the collision 
parameter, microwave resonant scans were simulated as described in Section 14.31 The 
simulated signal-to-noise ratio for T = 150 ns, 350 ns, 500 ns, 700 ns and 1000 ns was 
compared to experimental data measured at p = 250 mbar from [7] . It was assumed that 
T e = 2.51^ [20] and e\ = e<i = 55%. The upper and lower limits, determined from the 
1 sigma error associated with T,, are plotted alongside the previously measured data, 
shown in figure El 

At small T, Tj = 4.5 x 10 5 s _1 tends to underestimate the signal-to-noise ratio while 
Tj = 2.3 x 10 5 s" 1 overestimates the signal when T is large. At T = 200 ns both the 
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upper and lower limits are too small. This is due to the laser depopulation efficiency, 
which has a larger effect on the signal at short T. The experimental data displayed in 
figure El were measured during a different year (2006) to those data displayed in figured 
(2008). During this time the p flux was larger with higher stability and therefore E\ and 
€2 were likely to be higher than during 2008. 

6. Conclusions 

This study of the collision rates between pHe + atoms and the He medium has been used 
both to determine optimal conditions for microwave resonance profile measurements and 
to compare with theoretical predictions so that the system may be better understood. 
Experimental measurements of elastic collisions published in previous papers [19] 
prompted a re-evaluation of theory [20] which had originally overestimated the collision 
rate. Since the cross-sections of the two collision processes were predicted to be similar, a 
measurement of the inelastic collision rate was expected to introduce valuable knowledge 
about the interactions. 

For the first time inelastic spin exchange collisions have been been measured 
between the HF states of pHe + . Given the complexity of the system and the uncertainty 
in determining the initial parameters of the theoretical model, the measured values are 
in agreement with the recent theory [24] . More rigorous calculations are anticipated for 
a more thorough comparison. 

The laser depopulation efficiency is shown to be smaller than predicted but also to 
depend heavily on the conditions of the p beam and alignment. Simulations comparing 
the results to earlier data indicate agreement between the two methods but the inelastic 
collision rate is most likely to tend towards theory in the upper limits of experimental 
uncertainty. 

Collisional effects in El transitions have previously been shown to vary depending 
on the state measured, therefore other states are of interest. A study of the collision 
processes in the p 3 He + HF structure, which contains an additional degree of freedom 
due to the helion spin, is also planned for future work. 
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